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Abstract 

The  surface  regeneration  or  de-sulfurization  process  of  a  sulfur-poisoned  (i.e.  sulfur-covered)  nickel  surface  by  02  and  H20  has  been  studied 
using  first-principles  calculations  with  proper  thermodynamic  corrections.  While  02  is  more  effective  than  H20  in  removing  the  sulfur  atoms 
adsorbed  on  nickel  surface,  it  readily  reacts  with  the  regenerated  Ni  surface,  leading  to  over-oxidization  of  Ni.  Thus,  H20  appears  to  be  a  better 
choice  for  the  surface  regeneration  process.  In  reality,  however,  both  02  and  H20  may  be  present  under  fuel  cell  operating  conditions.  Accordingly, 
the  effects  of  the  partial  pressures  of  02  [po2]  and  H20  [pu2o]  as  well  as  the  ratio  of  po2/  pu2o  on  the  regeneration  of  a  sulfur-covered  Ni  surface 
without  over-oxidization  at  different  temperatures  are  systematically  examined  to  identify  the  best  conditions  for  regeneration  of  Ni-based  SOFC 
anodes  under  practical  conditions. 

©  2007  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

Sulfur-containing  compounds  are  the  most  commonly  encountered  contaminants  in  readily  available  fuels  and  are  difficult  to 
be  completely  removed  efficiently  and  economically.  Even  with  small  amount,  they  severely  degrade  the  performance  of  catalysts 
for  fuel  reforming  or  the  anodes  in  solid  oxide  fuel  cells  (SOFC).  The  poisoning  effect  or  the  degradation  in  performance  is 
generally  believed  due  to  the  adsorption  of  sulfur  on  the  surface  of  catalysts  or  SOFC  anodes,  which  blocks  the  active  sites 
for  desirable  reactions  [1].  Thus,  the  removal  of  the  atomically  adsorbed  sulfur  from  surfaces  is  a  critical  step  toward  sulfur 
tolerance. 

To  date,  the  de-sulfurization  of  various  materials  has  been  extensively  studied  via  experimental  and  computational  approaches, 
including  pure  metals  (Ni  [2-4],  Cu  [5],  Mo  [6,7]  and  Pt  [8,9])  as  electrodes,  small  clusters  (CoMo,  NiMo  [7,10],  N^P,  M02C 
[11],  M0S2  [12,13])  as  catalysts,  and  metal  oxides  (Ce02  [14-16],  Zr02  [17,18]  and  Fa203  [19])  as  electrolytes  in  fuel  cells.  It  is 
well  known  that  the  good  overlaps  between  sulfur  p  orbital  and  transition  metal  d  orbitals  result  in  strong  S-metal  bond  formation, 
making  the  de-sulfurization  process  difficult.  Small  clusters  of  several  metals  display  better  properties  for  de-sulfurization.  This 
enhancement  is  attributed  to  the  promoting  effect  in  CoMo  and  NiMo  cases  [7,10]  and  to  the  stabilization  of  nickel  d  orbital  in  the 
M2P  case  [11].  Among  the  metal  oxide  compounds,  a  highly  reactive  oxide  of  Fa2C>3  was  found  to  be  an  promising  material  to 
efficiently  remove  H2S  contaminant  from  fuels  due  to  the  dual  valance  charges  of  Fam  and  Fa1 * * IV  [19]. 

Furthermore,  various  chemical  reagents  have  been  thoroughly  studied  in  the  regeneration  process,  including  H2  [20-22],  O2 
[2,4, 16]  and  H2O  [15,23,24] .  The  surface  de-sulfurization  by  H2  forming  H2S,  which  is  the  reversed  reaction  of  the  highly  exothermic 
sulfurization  process,  H2S  -a-  S(a)  +  H2  [9,25-28],  has  shown  to  be  less  efficient.  On  the  other  hand,  the  oxidation  of  sulfur  by  O2 
and  H2O  to  form  SO2  shows  more  promising  results.  However,  over  dosing  these  reagents  on  metal  surfaces  results  in  metal 
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oxide  formation  and  degrades  the  catalytic  activity  of  Ni-based  catalysts.  Thus,  determining  the  right  conditions  to  de-sulfurize  a 
sulfur-poisoned  surface  without  over  dosing  is  a  critical  issue  to  achieve  surface  regeneration. 

Here  we  report  the  results  on  our  systematic  examination  of  proper  conditions  for  regeneration  of  sulfur-poisoned  Ni  surfaces.  In 
this  analysis,  we  used  DFT  to  compute  the  heat  of  reactions  and  thermodynamic  corrections  to  estimate  the  Gibbs  free  energies  of 
the  interfacial  reactions,  producing  phase  diagrams  important  to  the  de-sulfurization  processes  of  Ni  surface.  The  over-oxidization 
of  Ni  surfaces  by  excess  reagents  is  also  considered.  Further,  proper  conditions  for  the  use  of  O2,  H2O,  and  a  mixture  of  them  for 
de-sulfurization  without  over-oxidation  have  also  been  studied. 

2.  Computational  method 

The  adsorption  of  atomic  S(a)  on  Ni  surfaces  and  bulk  phases  of  pure  Ni  and  NiO  are  computed  by  Vienna  Ab  Initio  Simulation 
Package  (VASP)  [29-31],  implementing  the  density  functional  theory  (DFT).  The  exchange-correlation  function  treated  by  the 
generalized  gradient  approximation  with  the  Perdew-Burke-Ernzerhof  formulation,  known  as  GGA-PBE  [32,33],  is  applied  for  the 
total  energy  calculations.  Combining  the  accuracy  of  augmented  plane  waves  with  the  cost-effective  pseudopotentials  implemented 
in  VASP,  the  projector- augmented  wave  method  (PAW)  [34,35]  is  applied  in  the  basis  set.  The  kinetic  cut-off  energy  of  400  eV  is 
employed  to  ensure  the  convergence  with  respect  to  the  basis  set  and  the  efficiency  of  the  computational  cost.  The  Brillouin-zone 
(BZ)  integration  is  sampled  with  0.05  A-1  x  2  A-  1  spacing  in  reciprocal  space  by  Monkhorst-Pack  scheme  [36].  All  calculations 
are  performed  imposing  zero  magnetization. 

In  poly-crystalline  Ni  catalysts  or  electrodes,  Ni(l  00)  surface  forms  the  most  stable  S(a)  adsorption  and  is  considered  as  the 
most  vulnerable  surface  orientation  to  be  poisoned  under  SOFC  operation  conditions  [37].  Thus,  the  (1  0  0)  surface  is  applied  in  the 
current  calculations.  Although  sulfur  may  adsorb  on  surfaces  of  other  orientations,  it  is  easier  to  remove  sulfur  from  these  surfaces 
than  from  the  (100)  surface.  The  Ni  bulk  is  initially  computed  to  predict  the  theoretical  lattice  parameter.  The  result  shows  that 
the  computed  lattice  constant  of  3.58  A  is  in  good  agreement  with  the  experimental  value  of  3.52  A  [38],  and  employed  for  the 
surface  calculations.  The  surface  is  modeled  by  the  super  cells  of  four-layer  Ni(l  00)  slab  with  64  Ni  atoms,  as  shown  in  Fig.  1, 
followed  by  a  10- A  vacuum  space.  The  S(a)  adsorption  energies  have  been  initially  computed  to  test  the  convergence.  The  result 
show  that  the  surface  models  with  5  and  6  Ni(l  0  0)  slabs  or  larger  vacuum  space,  15  A,  have  similar  adsorption  energies,  less  than 
0.01-eV  difference,  to  the  applied  surface  model,  four-layer  Ni(l  00)  slab  and  10-A  vacuum  space.  The  Monkhorst-Pack  k-points 
are  set  as  2  x  2  x  1  along  (1  10)x(110)x(001)  directions.  The  modeled  surface  with  the  area  of  9.90  A  x  9.90  A  and  16  fourfold 
hollow  sites  (the  most  stable  adsorption  sites)  is  large  enough  to  avoid  any  adsorbate-adsorbate  interaction  and  to  simulate  minimum 
coverage  of  S(a)  with  the  highest  adsorption  energies.  The  top  three  layers  are  relaxed  and  the  bottom  one  layer  is  fixed  at  the 
computed  lattice  constant  to  represent  the  semi-infinite  bulk  crystal  beneath  the  surface  during  the  calculation. 


Fig.  1 .  The  Ni(l  0  0)  surface  is  constructed  by  64  Ni  atoms  with  a  10.4- A  vacuum  space.  The  blue  balls  are  the  surface  Ni  atoms.  The  dash  lines  represent  the  volume 
of  the  super  cell.  (For  interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  the  article.) 
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Considering  the  effect  of  temperature  and  partial  pressure  of  various  gases,  we  calculated  the  Gibbs  free  energy  G(T ,  p)  of  the 
whole  system  by  applying  the  thermodynamic  corrections  on  the  DFT  computed  energy,  £DFT,  which  represents  the  Gibbs  free 
energy  at  0  K  in  vacuum  or  the  Helmholtz  free  energy,  F(7),  at  0  K,  EDFT  =  G(0,  0)  =  F( 0). 

For  a  general  gas-surface  reaction: 

surface  1  +  gas  1(g)  — >►  surface2  +  gas2(g)  (1) 

where  surface  1  and  surface2  represent  the  reactant  and  product,  respectively,  in  the  condense  phase  whereas  gas  1(g)  and  gas2(g) 
correspond  to  the  reactant  and  product,  respectively,  in  the  gas  phase.  The  Gibbs  free  energy  for  reaction  (1)  can  be  expressed  as 


A G(T,  p)  —  (Ggas2  “h  Gsurface2)  (Ggasl  +  Gsurfacel)  (2) 

It  has  been  proven  [39,40]  that  the  Gibbs  free  energies  of  surfaces  have  relatively  small  variation,  <10meV,  in  a  wide  range 
of  temperature  (<1500  K)  and  pressure  (<100  atm)  and  the  vibrational  and  entropic  contributions  to  Gibbs  free  energies  of  atomic 
adsobates  are  negligible.  Therefore,  the  Gibbs  free  energies  of  solid  phases  can  be  approximated  by  the  internal  energy  computed 
from  DFT  calculations,  GsurfaCe(^  P)  ~  Gsurface(0,  0)  =  £^rfaCe-  contrasb  the  Gibbs  free  energies  of  the  gas  phases  are  strongly 
affected  by  the  temperatures  and  pressures  and  can  be  expressed  as 

Ggas(r,  p)  =  Ggas( 0,  0)  +  AGgas(0  T,  P°)  +  AGgas(7;  p°  -*  p)  =  E™  +  Hgas(T,  p°)  +  RT  In (Pgas/p°)  (3) 

The  first  term  is  directly  from  DFT  calculations.  The  second  term,  which  can  be  found  in  the  empirical  thermodynamic  database 
[41],  is  the  standard  enthalpies  of  the  gas  phase  molecules  contributing  from  rotations,  vibrations,  and  ideal-gas  entropy  at  p°  and  is 
temperature-dependent.  The  third  term  is  a  partial  pressure-dependent  function  from  Maxwell  relation  (where  R  is  the  gas  constant). 
Therefore,  at  a  given  pressure  and  temperature,  the  Gibbs  free  energy  can  be  obtained  by  adding  the  thermodynamic  corrections  in 
the  last  two  terms  to  the  DFT  results,  as  discussed  previously  [12,42,43],  and  the  Gibbs  free  energy  for  reaction  (1)  can  be  re-written 
as 


A G(T,  P)  =  (£?£ Le2  -  £?urfacel)  +  "  O  +  (#gas2  -  #gasl)  +  RT  (in  [pg as2J  -  In  |_PgaslJ  ) 


=  A£'dft(0,  0)  +  A H(T,  p°)  +  RT  In 


/7gas2 
_PgSLSl  _ 


(4) 


The  feasibility  of  the  gas-surface  reactions  can  be  determined  from  the  computed  A  G(T,  p ),  as  a  function  of  A£DFT,  enthalpy, 
temperature  and  pressure. 

The  thermodynamic  correction  of  external  electric  field  is  excluded  in  the  current  calculations  since  the  external  electric  field, 
which  exists  in  the  electrochemical  reactions  under  SOFC  operation  conditions,  has  relatively  smaller  contribution  in  the  Gibbs  free 
energy  calculations,  as  examined  in  the  Ni-S  system  [37]. 


3.  Results  and  discussion 


3.1.  Regeneration  by  O2  or  H2O 

The  sulfur  poisoning  effect  of  anode  for  an  SOFC  is  resulted  from  the  atomic  sulfur  adsorbed  on  the  Ni  surface,  which  blocks 
the  active  sites  for  fuel  oxidation  [1,37].  To  regenerate  the  surface,  therefore,  the  adsorbed  atomic  sulfur  must  be  removed  from  the 
surfaces  by  introducing  oxidizing  reagents  to  react  with  the  adspecies.  The  de-sulfurization  process  using  the  oxidant  reagents  of 
O2  and  H2O  can  be  described  as  follows: 

02(g)  +  S(a)  — >  SC>2(g)  +  Ni(l  00),  A£5DFr  =  -0.13eV  (5) 

2H20(g)  +  S(a)  -*  S02(g)  +  Ni(l  00)  +  2H2(g),  Afi?13'1'  =  4.66eV  (6) 

where  A£^FTand  A  Eg171  are  the  heat  of  reactions  (5)  and  (6),  respectively,  which  can  be  determined  from  DFT  calculations. 

Ah5  —  (£S02  +  ^Ni(100)l  _  h02  +  ^S(a) 

A^FT  =  (£dft  +  £dfto0)  +  2Euft)  _  (2£DFT  +  £DFT}  (8) 

As  described  in  the  S-Ni  system,  the  most  stable  S(a)  adsorbed  on  the  fourfold  hollow  site  of  Ni(l  0  0)  surface  with  the  lowest 
surface  coverage  [37]  considered  as  the  decisive  poison  species,  is  applied  in  the  present  calculation.  Since  the  non-adsorbed  species 
are  in  the  gas  phase  in  the  computed  temperature  range,  500-1600  K,  the  Gibbs  free  energies  for  reactions  (5)  and  (6)  are  given  by, 

Po2 


A G5  =  (Gso2  +  GNi(i00))  -  ( Go2  +  GS(a))  =  A£!?ft(o5  0)  +  A HS(T,  p°)  -  RT  In 


.P so2. 


(9) 
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A G6  =  (GSo2  +  Gm i  oo)  +  2Gr2)  -  (2Gh2o  +  Gm)  =  AE^iO,  0)  +  AH6(T,  p°)  -  2 RT  In 


PHjO 

Ph2Pso2 


(10) 


where  AHs(T,p°)  and  AH^,(T,p°)  are  the  standard  enthalpies  of  the  gas  phase  species  involved  in  reactions  (5)  and  (6),  respectively. 


A HS(T,  p°)  =  Hs02(T,  p°)  -  H02(T,  p°)  (11) 

A H6(T,  p°)  =  HSo2(T,  p°)  +  2Hn2(T,  p°)  -  2HH2o(T,  p°)  (12) 


By  computing  the  Gibbs  free  energies  of  AGs  and  AG 6,  as  functions  of  temperature  and  pressures  of  pso2 ,  Po2,  Pu2  and  pu2o , 
the  required  de-sulfurization  condition  can  be  expressed  in  the  thermodynamic  diagrams  in  Fig.  2.  The  yellow  regions  represent 
the  poisoned  phase  of  S(a)  adsorption  and  white  regions  represent  the  phase  of  clean  Ni  surface,  regenerated  by  O2  or  H2O.  The 
boundaries  between  the  yellow  and  the  white  regions  correspond  to  AG5  =  0  in  Eq.  (9)  and  AG 6  =  0  in  Eq.  (10). 


log 


PO2 
. Pso2 


=  0.434 


ae!?ft  +  ah5 

RT 


=  0.434 


( 77DFT  1  77DFT 

t^so2  +  E 


Ni(l  00)' 


)  -  (£Rft  +  £RFI  )  +  (HSo2  ~  Hq2) 


-02 


7DFT 


-S(a) 


RT 


(13) 


Ph2  O 

=  0.434 

'ae^  +  ahJ 

_Ph2Psq2_ 

2  RT 

=  0.434 

'(E°P  +  £DfJ00)  +  2E™f  )  -  (2Egf  +  £DFT)  +  +  1H^  _  2H^0 

2  RT 


For  the  O2  case  in  Fig.  2(a),  higher  temperature  requires  more  O2,  or  higher  log  [po2/ p$o2],  to  remove  S(a)  from  the  exothermic 
reaction  (5).  For  H2O  case  in  Fig.  2(b),  less  H2O,  or  lower  log  [pn2o/ Pn2Pso2],  is  needed  at  higher  temperatures  since  the  removal 
of  S(a)  by  H2O  in  reaction  (6)  has  a  high  endothermicity. 

It  is  noted  that  at  higher  or  lower  temperatures,  the  bulk  nickel  sulfide  will  be  formed  on  the  surface,  as  detailed  in  a  previous 
study  [37].  The  boundaries  between  S(a)  adsorption  and  nickel  sulfide  bulk  are  excluded  since  the  S(a)  adsorption  is  the  most  critical 
phase  in  the  poisoning  behavior;  thus  only  the  boundary  between  clean  Ni  surface  and  the  S(a)  adsorption  phase  is  considered  in  the 
this  analysis. 

Considering  the  pressure  effects  at  a  given  temperature,  higher  pressures  of  O2  or  H2O  can  clean  the  S(a)-poisoned  surface  more 
efficiently;  however,  exposure  to  excess  amount  of  O2  or  H2O  will  over-oxidize  the  regenerated  Ni,  resulting  in  O(a)  adsorption  and 
NiO  formation  [2,4].  While  O(a)  adsorption  can  be  readily  removed  by  H2  or  other  fuels,  the  formation  of  bulk  NiO  will  significantly 
de-activate  the  Ni  catalyst,  due  mostly  to  the  volume  change  associated  with  the  redox  reaction  of  Ni.  Taking  this  into  consideration, 
we  have  to  carefully  examine  the  Ni-NiO  bulk  phase  equilibrium  in  the  O2-N1  and  H2O-N1  systems. 

The  reactions  between  Ni  and  O2  or  H2O  forming  NiO  can  be  expressed  as 

\02(g)  +  Ni  -*  NiO,  A£°ft  =  -1.43  eV  (15) 

H20  +  Ni  NiO  +  H2,  A£?6FT  =  0.97 eV  (16) 


Fig.  2.  The  thermodynamic  diagram  for  S(a)  removal  process  by  (a)  02  and  (b)  H20  in  reactions  (5)  and  (6),  respectively.  The  yellow  and  white  regions  represent 
the  poisoned  S(a)  adsorption  and  regenerated  Ni  surfaces,  respectively.  (For  interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader  is  referred  to  the 
web  version  of  the  article.) 
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Fig.  3.  The  thermodynamic  diagram  for  (a)  O2  and  (b)  H2O  interacting  with  Ni  forming  NiO.  The  white  regions  correspond  to  the  phase  of  clean  Ni.  The  red  regions 
correspond  to  the  phase  of  NiO  bulk  formation  in  reactions  (15)  and  (16).  (For  interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader  is  referred  to 
the  web  version  of  the  article.) 


Similarly,  the  corresponding  Gibbs  free  energies  for  reactions  (15)  and  (16)  can  be  computed  as 
AGis  =  (£^fj  -  fT)  -  +  Ho2(T,  p°)  -  RTln[po2 ]) 


(17) 


AG 


16 


?DFT 

-NiO 


17J 

—  E 


-DFT 

Ni 


+  E™f  -  +  Hh2(T.  p")  -  HH2o(7;  p")  -  RT  In 


rDFT 


PH2Q 

Pn2 


(18) 


Shown  in  Fig.  3  are  the  thermodynamic  diagrams  for  the  02-Ni  and  F^O-Ni  systems  constructed  by  the  Gibbs  free  energy 
calculations.  The  red  regions  represent  the  over-oxidized  phase  of  NiO  formation  and  white  regions  represent  the  phase  of  pure  Ni. 
The  boundaries  between  the  red  and  the  white  regions  in  Fig.  3(a)  and  (b)  correspond  to  AG15  =  0  in  Eq.  (17),  and  A Gm  =  0  in  Eq. 
08). 


log[po2]  =  0-434 


n  77DFT  /  77DFT  1  9  77DFT \  tt 

ZjfcNiO  +  ZjfcNi  1 


RT 


09) 


log 


PWiO 
.  Pu2  . 


0.434 


7  DFT 
"NiO 


Z7DFT  1  77 

-  ANi  +  t 


DFT 

H2 


-  +  HHl  -  HHiq 


RT 


(20) 


Since  the  NiO  formation  from  O2  reacting  with  Ni  are  exothermic,  the  NiO  phases  appear  in  the  lower  temperature  range  at  a 
given  O2  partial  pressure  in  Fig.  3(a).  In  contrast,  the  NiO  formation  from  H2O  reacting  with  Ni  is  endothermic.  Accordingly,  the 
NiO  phases  appear  in  the  high  temperature  region  at  a  given  log[/7H2o/JPH2]  in  Fig-  2(b).  For  the  pressure  effect,  the  NiO  phases 
appear  in  the  region  with  high  O2  and  H2O  partial  pressures,  larger  values  of  log[po2 ]  and  log[/7H2o/JPH2]’  respectively,  at  a  given 
temperature. 

Therefore,  the  proper  condition  to  regenerate  the  surface  without  over-oxidization  can  be  determined  by  combining  Figs.  2  and  3, 
as  summarized  in  Fig.  4.  Insufficient  oxidant  reagents  of  O2  and  H2O  results  the  S(a)  phase  in  the  upper  left  and  lower  left  yellow 
regions  in  Fig.  4(a)  and  (b),  respectively.  The  over  dosed  reagents  of  O2  and  H2O  results  NiO  formation  in  the  lower  right  and  upper 
right  red  regions  in  Fig.  4(a)  and  (b),  respectively.  The  appropriate  conditions  to  regenerate  the  Ni  surface  are  located  in  the  center 
white  regions. 

Specifically,  the  partial  pressure  range  of  O2  in  the  white  (clean  Ni)  region  is  resulted  from  (19)  to  (13): 


A  log[po2]  =  0.434 


9/7DFT 

ZjGNiO 


-2E 


DFT 

Ni 


77 DFT  _  17DFT 
^S02  ^Ni(lOO) 


RT 


+  +  h*o2 


-  log[/?so2] 


(21) 


The  partial  pressure  range  of  [ph2o/th2]  in  the  white  region  is  resulted  from  (20)  to  (14): 


A  log 


Th2  o 
.  7*h2  . 


0.434 


9  77DFT  9  17DFT  77DFT 

^NiO  Zj^Ni  ^S02 


2  RT 


fdft 

^Ni(lOO) 


+  Hso2 


1 

21 °g[Pso2 ] 


(22) 


Therefore,  the  allowed  ranges  of  A  log[po2]  and  A  log[/?H2o//?H2]  for  the  O2  and  H2O  regenerations,  respectively,  depend  strongly 
on  the  partial  pressures  of  pso2-  the  lower  pso2  in  the  system  will  allow  wider  ranges  of  Alog[po2]  or  A  log[/?H2o/j7H2L  The 
partial  pressures  of  pso2  produced  from  the  O2  in  reaction  (5)  is  expected  to  be  higher  than  that  produced  from  H2O  in  reaction  (6) 
since  the  stronger  O2  oxidant  generates  SO2  faster  and  leave  more  residual  SO2.  Therefore,  the  last  terms  in  Eqs.  (21)  and  (22)  will 
further  narrow  the  proper  pressure  range  for  O2  but  widen  that  for  H2O  to  achieve  complete  regeneration  but  without  over  oxidation 
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Fig.  4.  The  complete  ( T ,  p)  diagram  for  (a)  O2  and  (b)  H2O  interacting  with  sulfur  adsorbed  Ni  surface.  The  yellow,  white,  and  red  regions  represent  the  phases  of 
S(a)  adsorption,  clean  Ni,  and  NiO  bulk,  respectively.  The  proper  conditions  for  surface  regeneration  without  oxide  formation  are  located  in  the  central  white  regions. 
The  boundaries  between  yellow  and  white  regions  correspond  to  the  axes  of  (a)  log[/7o2//7so2]  and  (b)  log[/?H2o/PH2Pso2b  boundaries  between  white  and  red 
regions  correspond  to  the  axes  of  (a)  log[po2]  and  (b)  log[/?H2o/pH2L  as  arrowed.  (For  interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader  is 
referred  to  the  web  version  of  the  article.) 


of  Ni.  In  addition,  the  pressure  range  of  O2  is  relatively  difficult  to  be  prepared  in  practical  experiments.  Thus,  H2O  is  considered 
as  a  better  reagent  for  regeneration  than  O2  as  observed  in  the  previous  experiments  [1,15,23]. 


3.2.  Regeneration  by  a  mixture  of  O 2  and  H2O 


Under  practical  fuel  cell  operation  conditions,  both  O2  and  H2O  may  be  present  in  the  system.  The  surface  de-sulfurization  and 
over-oxidization  can  be  described  as  follows: 


jc02  +  2(1  -  x)H20  +  S(a)  ->  S02  +  2(1  -  *)H2  +  Ni(l  00),  A E%fl  =  -0.13*  +  4.66(1  -  x)  (23) 

§02  +  (1  -  x)H2Q  +  Ni  ^  NiO  +  (1  -  x)H2,  A£g8FT  =  -1.43x  +  0.97(1  -  x)  (24) 


where  *  varies  from  0  to  1 ,  depending  on  the  ratio  between  O2  and  H2O.  At  *  =  0,  it  is  the  pure  H2O  regeneration  process  as  described 
by  reactions  (6)  and  (16)  while  at  *=  1,  it  is  the  pure  O2  regeneration  process  in  reactions  (5)  and  (15). 

The  Gibbs  free  energies  of  reactions  (23)  and  (24)  can  be  expressed  as 


A  G23  =  A£'?3ft  +  A  H23  +  RT  In 


2(1 -x) 

pso2Pk2 

x  2(1  -x) 

Po2Ph2o 


A G24  =  AEW1  +  A H24  +  RT  In 


(1 — JC) 

PUi 

r/2  (1 — v) 

Po;Ph2o 


(25) 

(26) 


Where  the  internal  energies  and  enthalpies  are 

A E™1  =  2(1  -  x)E™f  +  E™  +  E^m)  -  +  2(1  -  x)E™  +  E™$) 

AH23  =  2(1  -  x)HH2  +  HSo2  -  (xHq2  +  2(1  -  x)HH2o) 

a  77DFT  / 1  \  77DFT  1  rDFT  ( x  77DFT  1/1  77DFT  1  rDFT^ 

Ah24  -  (1  -  X)t h2  +  tmo  -  y2t02  +  I1  “  x)^H20  +  ^Ni  ) 

A H24  =  (1  -  x)Hu2  -  (f  Ho2  +  (1  -  x)Hu2o) 


(27) 

(28) 

(29) 

(30) 


From  the  Gibbs  free  energy  calculation,  the  boundaries  among  the  phases  of  clean  Ni,  S(a)  and  NiO  from  different  values  of  *  is 
shown  in  Fig.  5. 

x/2  (T—  x)  0  s  n  —  jc) 

For  the  pressure  effect,  the  required  reagents  to  remove  S(a)  in  reaction  (23),  log[/?02  /?h2Q  / ^so2^h2  ]’  an<^  t0  ^orm  NiO  in 

reaction  (24),  \og[pXQ^ decrease  as  the  v  increases  because  less  oxidants  are  needed  in  both  cases  as  the  component 
of  stronger  O2  oxidant  raises.  For  the  temperature  effect,  reactions  (23)  and  (24)  are  changed  from  exothermic  to  endothermic  as  * 
decreases  since  the  pure  O2  regeneration  (x=  1)  in  reactions  (5)  and  (15)  are  exothermic  and  the  pure  H2O  regeneration  (*  =  0)  in 
reactions  (6)  and  (16)  are  endothermic.  This  result  explains  the  boundary  change,  from  left  to  right  concave,  as  *  decreases. 
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Fig.  5.  The  ( T ,  p )  diagrams  for  (a)  S(a)  removal  and  (b)  NiO  formation  in  a  mixture  of  O2  and  H2O  of  different  ratio.  The  white,  yellow,  and  red  regions  correspond 
to  the  phases  of  clean  Ni,  S(a)  adsorption,  and  NiO  bulk,  respectively.  The  related  Gibbs  free  energy  calculations  for  (a)  and  (b)  are  listed  in  Eqs.  (25)  and  (26), 
respectively.  (For  interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  the  article.) 


In  addition,  shown  in  Fig.  6(a)  and  (b)  are  the  contours  for  the  required  reagent  pressures  to  de-sulfurize  the  surface  and  to  oxidize 
Ni,  respectively,  in  different  x  values.  The  proper  pressure  ranges  are  in  between  the  orders  of  lo g[Po2PH^/ ^so2^h2  ^  f°r  S(a) 

removing  and  1  oglPo^PH^/ Ph2  ^1  f°r  NiO  formation  (between  the  two  contours).  Following  the  analogous  works  in  Eqs.  (21) 
and  (22),  the  orders  of  proper  pressure  ranges  can  be  expressed  as 


A  log 
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(31) 


The  order  of  proper  pressure  range  is  independent  of  x  under  thermo-equilibrium  conditions.  It  means  that,  in  the  two  extreme 
cases,  A  log[/?02  ]  (at  v=  1)  equals  to  A  log[pu2o ]  (at  x  =  0),  which  is  consistent  with  the  results  in  Eqs.  (21)  and  (22).  In  addition, 
considering  the  last  term  of  l/21og[/?so2]  in  Eq.  (31),  we  can  see  that  stronger  O2  oxidant  in  the  reaction  (23)  will  generate  SO2 
faster  and  leave  more  residual  SO2.  Therefore,  proper  range  for  A 1  °gVPo2P\{2o  / Pu2  ]  exPecte(i  t0  be  narrower  at  higher  x 
values,  which  is  also  consistent  with  the  conclusion  in  Section  3.1. 

It  is  noted  that  our  analysis  is  applicable  only  to  thermodynamic  equilibrium  conditions,  predicting  the  ultimate  equilibrium  states 
under  a  given  set  of  conditions  (temperature  and  pressure).  The  kinetics  of  these  reactions  are  important  as  well.  However,  detailed 
kinetic  study  requires  additional  calculations  of  the  intermediate  and  transition  states  in  the  reaction  paths  (potential  energy  surface) 
and  molecular  dynamic  simulations.  In  addition,  surface  morphology  might  change  due  to  nickel  atoms  sintering  or  migration  on  the 
surface  during  the  regeneration  process.  Further,  oxygen  ions  transported  from  the  cathode  side  may  interact  with  the  adsorbed  sulfur 
on  Ni  surfaces  under  practical  fuel  cell  operating  conditions.  These  dynamic  effects,  which  have  not  been  taken  into  consideration  in 
this  study,  may  significantly  influence  the  regeneration  process.  The  methodologies  to  address  these  issues  are  still  under  investigation 
and  will  be  reported  in  subsequent  communications. 


Fig.  6.  Contour  plots  of  the  required  pressure  ratios  for  (a)  log [ p'q^ P(\\2o  I ^so2  ^h2  m  S(a)  removal  reaction  and  (b)  1  og [ px(!^ p^ q  Vph2  A)]  in  NiO  formation. 
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4.  Conclusion 

Regeneration  of  sulfur-poisoned  Ni  surfaces  by  O2,  H2O,  and  a  mixture  of  them  has  been  examined  using  ab  initio  atomistic 
thermodynamics.  In  the  regeneration,  both  the  removal  of  the  atomically  adsorbed  sulfur  and  over-oxidation  of  Ni  surfaces  to  form 
NiO  are  considered.  The  required  temperatures  and  pressures  of  the  chemical  reagents  needed  to  remove  the  atomic  sulfur  adsorbate 
and  to  from  NiO  are  predicted  and  presented  in  thermodynamic  diagrams.  While  both  O2  and  H2O  may  be  used  for  the  regeneration 
process,  H2O  appears  to  be  a  better  option  because  of  its  broader  pressure  range  for  removing  surface  sulfur  without  over-oxidizing 
the  Ni  surface. 
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